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ABSTRACT 

The problem of a periodic scalar field on a two-dimensional dynamical random 
lattice is studied with the inclusion of vortices in the action. Using a random 
matrix formulation, in the continuum limit for genus zero surfaces the partition 
function is found exactly, as a function of the chemical potential for vortices of unit 
winding number, at a specific radius in the plasma phase. This solution is used 
to describe the Kosterlitz-Thouless phenomenon in the presence of 2D quantum 
gravity as one passes from the ultra-violet to the infra-red. 
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1. Introduction 



The Kosterlitz-Thouless (KT) phenomenon [1, 2, 3] has had far reaching implications for 
our understanding of effectively two-dimensional (2D) physical systems. Its myriad applica- 
tions include those to 2D superfluids [1], crystal growth [4], lattice gauge theories [5], and 
superstring theories [6]. While this phenomenon is usually studied on a fixed regular lattice, 
it is known that many aspects of 2D statistical mechanics become more tractable when stud- 
ied on lattices for which local connectivity is itself a dynamical variable, principally because 
they may be mapped to exactly solvable problems of random matrices [7] . In the language of 
two-dimensional quantum field theory this corresponds to the introduction of 2D quantum 
gravity. They seem to display the same critical phenomena as on fixed regular lattices but 
with mild modification of critical exponents; for example, the Ising model transition is third 
order rather than second [8]. The possibility of finding some non-perturabtive results for the 
KT vortex plasma motivates a study of this phenomenon on such fluid random surfaces. In 
this paper, using the example of a circular scalar field at a particular radius, it is shown by 
an exact solution in the continuum limit how the gas of unbound vortices leads to topologi- 
cal order in two dimensions, making no assumptions about diluteness save the restriction to 
vortex winding number one (the most relevant). 

The random matrix model employed is a particular case of those studied in [11] which may 
be reformulated as an 0(2) model on a dynamical triangulation [18, 20]. On a single oriented 
link with periodic boundary conditions, representing the circular target space, one has an 
NxN complex matrix M — the models are the non-critical string analogue of Weingarten's 
matrix model [12] — the full partition function being 



Z= DM exp ( - Tr[M^M] + — Tr[M^MM^M + MMM^M^} + h^NTr[M + M f ] 



Expansion in k, the bare cosmological constant, generates a dynamical quadrangulation 
(square simplices glued pairwise along edges) with the edges embedded on the periodic link; 
also included in the action is a coupling h to vortices (and anti-vortices) of unit winding 
number, which cut small holes in the random surface that wind around the target circle. 
Gross and Klebanov have shown [15] that, in the double scaling limit [14], the vortex-free 
c = 1 matrix model partition function defined on a lattice is equivalent to that on the 
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continuous real line provided the lattice spacing a is less than some known critical value a c . 
At phase transition occurs (believed to be of KT type) due to the appearance of 

a marginal operator. The same holds true of the complex matrix models of ref. [11], from 
which one may deduce that (1) represents the case a = a c ; this means that the periodic link 
has radius r = a c /2ir = |r s d, where r s d is the self-dual radius. This places one deep in the 
purported vortex-plasma phase since unit charge vortices become relevant for r < 2r s ^. 

It is precisely at r = \r s d that one may reformulate the problem in terms of a solvable 
0(2) model by introducing an inducing field [20, 18]. Using an NxN Hermitian matrix 
variable one can rewrite (1) as 

Z oc J DM D<p exp ^- Tr[M f M] - Tr 2 + yj^Tr[M^M<j> + MM f 0] + hVNTr[M + M f ] j , 

(2) 

which after diagonalising and performing the Gaussian M-integrals gives 

r N 

z(k 6o) = / n dX * ?m-nJ2 2A ' + J2 lo § ( A * - a j) - J2 lo § ( 26 ° - x * - x ^ 

»=1 * i+3 id (3) 

+ ANh 2 J2(bo -A l )- 1 ) . 



In eq.(3) A are the eigenvalues of 0, and there has been a redefinition v2/t = l/4&o and 
rescaling A — > 2\\fN and /i 2 — > 8h 2 /bQ. The case h = has been studied extensively by 
I.Kostov and collaborators [19] as a loop gas problem, who found c = 1 critical exponents. 
Thus despite the presence of a marginal operator, at the radius chosen, over which one has 
little control, one can still study the KT phenomenon as h is turned on since the ultra-violet 
(UV) limit, corresponding to h = as we shall see shortly, is that of a massless scalar theory. 

At h 0, the vortex-anti- vortex pairs appear as an extra singular term in the potential 
for this eigenvalue problem, which at large iV may be treated by saddle-point methods [21]. 
Indeed, the operator coupling to h was identified already in the loop gas approach [19] and 
has a clear geometrical meaning. If one expands Z (2) in Feynman diagrams, a typical piece 
of planar diagram will appear as in Fig.l. A "bug" crawling on the surface goes 2n around 
the circle every time it crosses an arrowed propagator < M^M >. Closed loops of these (the 
loop gas) delineate regions of constant scalar field, X say, which since the target lattice has 
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only one link, takes only one value. Vortex-anti-vortex pairs appear connected by cuts in 
the X field. If the length of a cut is /, then in terms of the field it acts like a hole of length 
21 on the surface with one marked point on the boundary. The other marked point, given 
by the other vortex, is not independent since the two are constrained to be length / apart. 
This is the origin of the macroscopic loop operator coupling to h 2 in (3) . 

Thus the Boltzmann weights on a globally connected planar lattice take the form 

logZ(6 , h) = h* vortice % tot&l lensth of loops & cuts . (4) 

graphs 

One may therefore think of the arrowed lines as flux tubes which, because of the peculiar 
single valuedness of the target lattice, carry a unit charge's worth of electric flux for a 
Coulomb field. — log h is the chemical potential for vortices, with small h corresponding to 
the dilute regime. Flux tubes densely populate the surface and the continuum limit for area 
is achieved when they become infinite in length as bo is tuned to its critical value [19]. 

2. Saddle-Point Solution. 



Let us now study the effect of turning on h by finding the continuum limit of (3) . The 
leading order of the large N limit of Z(h, bo) describes surfaces of spherical topology. There 
exists a systematic approach to computing the 1/N corrections [18], corresponding to higher 
genus surfaces, but here we will concentrate on genus zero. Introducing an eigenvalue density 
p(\) for the problem (3) , the saddle-point equation for p is 

| d , pW (p^ + _^). 2A __^l p , (5) 

a 

where the support of p is on the interval [a, b}. This integral equation has been solved for 
h = by Gaudin [22]. The same method will be used to solve h ^ also. To transform (5) 
to an integral equation with Cauchy kernel make the following redefinitions 

A -> bo - VA + BX , A + B = (b - a) 2 , A - B = (b - b) 2 . (6) 
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Then 



/i oh 2 
dli P(») p x^ = 2 (^+^ - bo) + ^-^ • (7) 

-l 

The inversion formula for this problem is [23] 

A and B are determined by the normalisation and positivity conditions for p; 

l 



ft2 

V^i + 5^ - 6 + 



(8) 



l 



(9) 



To study the critical behaviour one need only know the form of p (8) for A — > — 1, though 
it may actually be found completely in terms of elliptic integrals. The continuum limit of 
surfaces is achieved as the bare cosmological constant is tuned to its critical value, given by 
the condition b = bo, when the end of the support of p meets the singularities in (5) . After 
some algebra, and performing the h-term integral, one can rewrite the solution (8) in the 
form 

P(A) = (BMW - (A + BA)n(A) - 2{a+ IX^b ) > < 10 > 

where 

n(A) = 7—J* = p ( - — - — ) , (ii) 

J Vl-^sin 2 ^ VA-cos20; 

Q V 

an elliptic integral of the third kind, while k 2 = 2B/(A + B), p = cos 20, and K (E) denotes 
the complete elliptic integral of the 1st (2nd) kind. From the positivity condition one finds 

nBh : 



2E(k)VATB - vrfeo + ^ ^ _ ^ = . (12) 
Using (10) , the normalisation condition requires a little more manipulation before it can be 
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brought to the form 

9 2h 2 
-^(A + B)((l - k 2 )K 2 (k) - E\k)) - j^= {K{k) + n.s.) = 1 . (13) 
vr z ix\J A — B 

n.s. stands for terms non-singular in the limit b — > bo (k — > 1), which will be unimportant 
for the continuum limit. To approach it one introduces the cut-off 5 — > and renormalised 
parameters M and £ as follows [19]: 

— - 1 = M 2 S 2 

- (14) 



M _1 is the renormalised dieletric constant of the Coulomb gas, M the renormalised boundary 
cosmological constant coupling to the length of cuts in X. Returning to (10) and evaluating 
it perturbatively in 8 one eventually finds (£ > M) 

_ „^ IQy/A... .... /-s —j /i 2 v/£ 2 - M 2 . . 

M) = — ^(log M5) Ve 2 - M 2 - ^ V ^ 2M + • • • (15) 

where . . . stands for terms subleading as 5 — > 0. The bare vortex parameter ft, couples to 
a relevant operator and therefore must be tuned to zero for finite renormalised parameter. 
The scaling law follows from (15) by requiring the 2nd term to be of the same order as the 
first 

h 2 = -h 2 R 5 3 \ogM5 . (16) 

This agrees with KPZ scaling [16] at r = ^r s d but there is a logarithmic scaling violation; 
such violations are well-known at c = 1, occurring also for the cosmological constant [17]. 
The latter is determined in terms of M by expanding (12) and (13) in 5; 



7l7? 2 

7vb = 2v / 2A(l - M 2 5 2 log MS) + lo S MS + ■ ■ ■ 

1 = ^(1- 2M 2 5 2 log 2 MS) - ^^5 2 log 2 M5 + ... 
ix l AirM 

which, omitting + . . . now, determines A = 7r 2 /4 and 



(17) 



A* = (M 2 + ^M)log 2 M5, (18) 



where bo = y/2 + 5 2 Ar defines the renormalised cosmological constant conjugate to renor- 
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malised area of the surface Ar. 5 is thus the physical cutoff. 

One can now eliminate M in p (f 5) in terms of the quantities Ar, Hr. Apart from the 
logarithms, (18) is cubic in M and it is convenient to use the parametric solution 

6 



where 



cos# 



M = 2 



hi 



A 7T 

cos 

3 3 
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3/2 



A 



(19) 



(20) 



v^vr 2 ^Ay ' " log 2 M5 ' 
The equation (18) is sketched in Fig. 2. (19) corresponds to the largest positive root, this 
being the physical one such that M — > oo as Ar — > oo. (15) and (19) constitute the solution 
for the scaling part of p, from which physical quantities may be calculated in the continuum 
limit. 

The first thing to notice is that (18) , which is the positivity condition for p, has no 
solution for h\ > V2tt 2 (A/3) 3 / 2 . In this case the singular term in the potential (5) has 
become too strong and there is no smooth large N limit. Vortices contribute vacuum energy 
which can be negative ~ q^r a r an d so the partition function at fixed cosmological constant 
does not converge in area Ar if Ar < (Hr too large). 

As an intermediate step to calculating the partition function of spherical surfaces, it is 
convenient to first derive the macroscopic loop expectation at genus zero i.e. the partition 
function for surfaces with the disc topology. 



u>(z) = J dl Ld(l)e~ 



zl 



(21) 



where cu(l) is the renormalised partition function for discs of renormalised perimeter length 
I with one marked point on the boundary. Using (14) and (15) 



u(z) 



M 



Pit, M) 



= -\ogM5 ( ^z 2 - M 2 ( log 



+ Vz 2 - M 2 
M 



h 2 

V2- ' ;/ '' 



7T 2 MZ 2 



+ 



\h\ ( it I 



71 \2z 2 zM / 
(22) 

In deriving (22) individual terms in the leading order which are analytic in z have been 
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omitted since they correspond to surfaces of infinitesimal area bounded by finite perimeter 
[24]. (22) for Hr = was given in [19]. As a useful check one may verify that Co correctly 
solves the finite difference equation derived by Kostov directly from the continuum limit of 
(5) (see e.g. ref.[18]). The partition function for spherical surfaces may be derived from 
Q(0) = dhgZ/dh 2 , 



(0) = 2V2(logM5) (m+J&^ 

n l2A R ( n-0 1 cos# 
= 2 \ / cos 



(23) 



2 cos 2 ^ 



by integrating w.r.t. h 2 at constant feo- The constant of integration is c ~ A 2 ^ by dimensions, 
up to logarithms. Using (20) it is simpler to integrate w.r.t. 0, dh 2 = —(A/3) 3 ^ 2 V2ir 2 sin 0d0+ 



, „ Ai / e-ir 1, 0-7i 1 . A0-7T . 



This represents the continuum limit where all subleading terms + . . . have been dropped. 
Equations (22) and (24) are the main results of this paper, representing the continuum limit 
on genus zero surfaces for the vortex gas with arbitrary coupling. 

One can use these results to investigate the number of effective degrees of freedom in 
various regimes (see refs.[10, 9] for a similar non-perturbative study of the Sine-Gordan 
model coupled to 2D gravity). In particular one can use (24) to find the string susceptibility 
7str, given by the fixed area partition function log Z(Ar) ~ A^ 3+7str exp(A^yl^). The UV 
limit Ar — > occurs for A# — > oo, in which case (24) yields 7 s t r = 0, as appropriate for 
a massless scalar field. As the scale Ar is increased above that set by h R the screening 
effects of the vortex condensate come into play and the behaviour changes. To determine 
7str more generally one must identify the effective cosmo logical constant, given by the sum 
of contributions from gravity (A#) and matter (— A^), which is conjugate to area. This 
constant is zero at the critical point of (18) allowing identification of A^. Since C = one 
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finds 



M 2 c \og 2 5M c = ^ 



V(A^) 2 



log 3 5M c = l (25) 



h 2 ^ 1/3 

M r .= ' 



As usual, the presence of logarithms complicates A^ which moreover satifies a transcendental 
equation. A C R = 3M 2 neglecting logarithms. 

To find 7 s tr in the infra-red (IR) limit Ar — > oo one can proceed analytically. Defining 
Ar = A C R + Ar and M = M c + M, for small 9 one may expand all formulae for Ar <C 
A C R , M <C M c and find expressions for the small variables; even keeping track of logarithms 
is not too tedious. Expanding about (25) one finds 



log M c 5 y 6 y A^ 3 V 3 
U> 8* _ 3/3 \ J 1/lx . 1(i> „ . 2.1/ 



(26) 



— r + — - ... = -- — m; Ap log 2 5M C + + . . . 

2 4! 2 V A i?/ M c logM c 5 



and expanding Ar in 9 gives 



= 3^ + O(Af). (27) 



Expanding the partition function (24) in 6? one finds 

logZ= { f£ ( a + pe 2 + ^ + ... + (e + vA R + r,A 2 R + ...j) , (28) 
log d A A c R v / 

where a, /3, 7, . . . are known constants and the constant of integration c = e + . . . is analytic 
as Ar — > since it does not depend upon h, only Ar. Thus 

logZ = analytic + ^g- A * + 0(Af ) , (29) 
log 5 2 A c Ry /A c R 

confirming that 7 s t r = —1/2 in the IR limit, where the cuts in X disorder it to leave pure 
gravity. Note that the analytic terms as Ar — > come from the lower cutoff on the area 
integral, most importantly they do not contain logarithms of Ar. 
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The change in behaviour from the UV to the IR can be understood from the Coulomb 
gas in terms of screening effects [1]. At small areas A# — > oo and the renormalised M _1 — > 
(18) , all flux lines being of the order of the UV cutoff on the surface; this is the c = 1 regime 
where vortices are strongly bound to anti-vortices. At large areas A# — > A C R and M — > M c 
(25) . When Hr = 0, since the surface is dense with flux loops, to achieve large area the 
dielectric constant is tuned to oo (M c = 0). In the presence of the charge gas (h,R ^ 0) 
the inverse permittivity is M — M c , the (non-universal) shift M c being due to dielectric 
polarization. It clearly required the non-perturbative treatment to see this since (25) is not 
analytic as fiR — > 0. 



3. Conclusions. 



Unfortunately it appears difficult to conceive of a general method for obtaining exact 
results for r ^ \^sd or higher vortex winding number. However one may add the winding 
number ±2 operator, m2Tr[MM + M^M^], to (1) and the integrals are still Gaussian, 
yielding 



Z = J Y[d\i exp(J2 log (A, - Xj) - N ^ 2A? - \ log ((2&o - A, - A,; 



\2 / ~ \2 



far) 



i+i ' (30) 



+ 8iVmi^(6 -^ 2 -A i )- 1 ) . 



instead of (3) , with h = mi. The geometrical meaning of the last two terms in the action 
(30) is illustrated in terms of flux lines in Fig. 3. The saddle-point equation to be solved is 

b 

, / x fr, 1 1 1 1 1 N. ^ 2(mi) 2 

du p(u) P- + — : + — = 2A 

p pvp; y x-fi 2 2& - rh 2 - A - /i 2 2b + m 2 - A - fi J (bo - \vrt2 - A) 



2 



a 



(31) 

but this is more difficult than (5) ; the author has not managed to solve it to the extent 
that useful information can be gleaned, for example on questions of multicriticality as fri2 
is tuned as a function of rh\. One may however give a hand- wave argument that the IR 
limit is pure gravity. Consider mi = and m<i ^ in (31) for simplicity. According to the 
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discussion of the previous section, the scaling regime is described by (A' + 2 m 2)^ — bo — A 
and 

oo 

[ dn'pU') ( Ptt^ - — T-r— — - - 5-^ 1 = 2(SX' + 5-m$ - b ) , (32) 

J V ' V A' - p! 2 A' + p> 2 2mf + A' + p' J V 2 2 u; ' K J 

b> 

where m R = 1712/6 is consistent with KPZ scaling. The IR limit is equivalent to — > 00 
in which case the 3rd term in the kernel of (32) may be dropped, leaving the saddle point 
problem of the 0(1) model in Rostov's classification of the dense phase of 0(n) models [19] 
i.e., pure gravity. 

For the more general models of ref. [11] one can ask about the extent to which similar 
methods of solution may be applied, in particular the introduction of induction fields to 
obtain plaquettes from an action quadratic in M. A particularly interesting problem concerns 
the effects of an extrinsic curvature dependence, which can smooth badly behaved surfaces 
[13]. This may be included in the Weingarten-type [12] models by introducing new flavours 
of link matrix. For example at c = 2, two flavours M, N suffice since the worldsheet normals 
are ±1. The extrinsic curvature, given by ^(1 — cos#) = 0, 1 for parallel (9 = it) or anti- 
parallel (9 = 0) neighboring plaquettes where 9 is the angle subtended, can be included by 
using an action 

Tr[Mt(/)M(/) + N*(l)N(l)] + K Tr[iV t (/)M(/) + M^(l)N(l)} - k ^ Tr[P M + P N ] , (33) 

l plaq 

where J2i * s the sum over all links of a square lattice; Pm is the product of M's around 
plaquettes of orientation + and Pjv the product of iV's around plaquettes of orientation 
— ; K is the bare extrinsic curvature coupling. The c-dimensional generalisation has Z c 
symmetry and c(c — 1) flavours, and one assumes that the continuous symmetry would be 
restored at any critical point. These models are difficult to solve in full generality, but 
simplified versions which still incorporate extrinsic curvature can be reduced to eigenvalue 
problems. This will be dealt with in a future publication [25] 

To summarize, in this paper the matrix model of random surfaces introduced in ref. [11] 
has been further investigated and applied to the two-dimensional vortex plasma in the pres- 
ence of gravity. The exact solution on genus zero was found at half the self-dual radius as 
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a function of renormalised cosmological constant and vortex chemical potential (22) (24) . 
In general these matrix models seem to give new possibilities for exact solution where the 
traditional Hermitian models are too difficult, the principle improvements being: 

1. Vortex operators are manifest and easily manipulated. 

2. Exact solutions are possible even in the presence of non-tree-like embeddings. 

3. Extrinsic curvature can be trivially introduced and, in certain cases [25], leads also to 
solvable models. 

4. The matrix models have gauge symmetry [11] and are in fact closely related to conven- 
tional non-linear lattice gauge theories, for which there is an abundance of expertise 
available. 

The accompanying geometrical picture confirms the KT hypothesis of unbound vortices, 
which in the renormalised scheme are joined by Coulomb flux tubes, cuts in the periodic 
scalar field that disorder it at large length scales. Hopefully this approach can be further 
developed to extend the remarkable non-perturbative answers which matrix models furnish. 

Acknowledgements: I have benefitted from discussions with I.Klebanov, I.Kostov, M.Martin- 
Delgado, and M.Staudacher. Financial support was provided by S.E.R.C.(U.K.) post-doctoral 
fellowship RFO/B/91/9033. 

Note Added: The reader's attention is also drawn to recent numerical studies of the XY 
model on random surfaces [26, 27]. 
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FIGURE CAPTIONS 



Fig.l. A typical piece of planar Feynman graph for action (1) . A single vortex-anti- vortex 
pair and three "vacuum" flux loops are shown arrowed. Dotted line is < 00 > pro- 
pogator. 

Fig. 2. A sketch of the cubic equation (18) . 

Fig. 3. Vortex flux configurations corresponding to the expansion of action (30) in powers of 
mi and m<i- 
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